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Abstract 

We present a new semiclassical method that yields an approximation to the quantum mechanical 
wavefunction at a fixed, predetermined position. In the approach, a hierarchy of ODEs are solved 
along a trajectory with zero velocity. The new approximation is local, both literally and from a 
quantum mechanical point of view, in the sense that neighboring trajectories do not communicate 
with each other. The approach is readily extended to imaginary time propagation and is particu- 
larly useful for the calculation of quantities where only local information is required. We present 
two applications: the calculation of tunneling probabilities and the calculation of low energy eigen- 
values. In both applications we obtain excellent agrement with the exact quantum mechanics, with 
a single trajectory propagation. 
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I. INTRODUCTION 

The first semiclassical method, the WKB method|l|, was pubhshed almost simultaneously 
with the publication of the Schrodinger equation in 1926. Since then, semiclassical methods 
have continued to attract great interest for two primary reasons. First, semiclassical methods 
give insight into classical-quantum correspondence. Second, for large systems they hold 
the promise of significant computational advantages relative to full quantum mechanical 
calculations. In particular, in recent years much progress has been made in the chemical 
physics community in developing time- dependent semiclassical methods that are accurate 
and efficient for multidimensional systems. By semiclassical methods, one generally means 
the calculation of a quantum mechanical wavefunction or propagator via the propagation 
of classical (or classical-like) trajectories. Mathematically speaking, semiclassical methods 
cast the time-dependent Schrodinger equation (TDSE), which is a PDE, in terms of ODEs 
related to classical equations of motion. From a physical point of view, the semiclassical 
methods try to circumvent the non-locality of quantum mechanics. 

In this paper we present a new method that is to some extent a cross between a numerical 
grid method and a semiclassical method. The method is derived by inserting a trial form e*'^/^ 
into the TDSE and calculating the time-dependent complex phase at a stationary position. 
Since the method yields an approximate solution of the TDSE at a fixed position it is a cousin 
of grid methods. However, we compute S by solving a hierarchy of spatial derivatives of S 
along a stationary trajectory; information from the neighboring trajectories is incorporated 
only through the initial conditions. This property makes the approximation local from the 
quantum mechanical point of view, and hence a cousin of semiclassical methods. We refer 
to the new method as the zero-velocity complex action method (ZEVCA) since it employs 
a complex phase (action) and stationary trajectories. 

The substitution of e*'^/'^ as an ansatz in the TDSE is the same starting point as sev- 
eral methods that are based on the hydrodynamic and Bohmian formulations of quan- 
tum mechanics . More specifically, ZEVCA is related to the Derivative Propaga- 
tion Method j^, the Trajectory Stability MethodQ] and Bohmian Mechanics with Complex 
Action 81 . As in the case with ZEVCA, these methods incorporate equations of motion for 
a hierarchy of spatial derivatives of the phase (and amplitude) that are calculated along 
trajectories. But unlike ZEVCA, in these other approaches the time dependent trajectories 
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propagate in either real or complex configuration space. The ZEVCA formulation employs 
fixed trajectories that yield the time dependence of the wavefunction at a fixed, predeter- 
mined position in configuration space. In reference jsl section 7.2, Wyatt considers the 
solution of the global hydrodynamic equations of quantum mechanics on fixed grid points 
(Eulerian grid) but he dismisses its usefulness as a numerical tool. The ZEVCA formulation 
shows how to obtain useful output of a local propagation at a single grid point. 

For a number of quantum mechanical calculations, such as thermal rates or tunneling 
probabilities, knowledge of the wavefunction in all of configuration space is unnecessary: 
these quantities can be calculated by data at a single position or a small interval of space. 
For such calculations, ZEVCA has a significant numerical advantage since it produces local 
information at a predetermined position. The first application of ZEVCA that we present in 
this paper is the calculation of tunneling probabilities; ZEVCA is well suited for this calcu- 
lation since tunneling probabilities can be calculated from a time integral of the probability 
current at a fixed position. The second application is the calculation of the lowest energy 
eigenvalue of a bound potential. For this application we make a minor modification to the 
ZEVCA formulation to adapt it to imaginary time propagation. In both applications we 
calculate the results by propagating just a single trajectory. 

This paper is organized as follows. In section [TTl we formulate ZEVCA for the real time 
Schrodinger propagator (section III Al) and the imaginary time Schrodinger propagator (sec- 
tion [HE])- Since the multi-dimensional representation of ZEVCA is somewhat more com- 
plicated, we present only the one- dimensional derivation. Section [III] is dedicated to the 
implementation of ZEVCA for the derivation of tunneling probabilities (IIII Al) and the first 
energy eigenvalue of a bound potential (jlll B|) . In section HVl we present our summary and 



concluding remarks. 



II. FORMULATION 



A. The ZEVCA real time Schrodinger propagator 

We start by inserting the ansatz [ol, [lo|, 11| 



ipi^x, t) = exp 



^S{x,t) 



(2.1) 
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into the TDSE 

thtfjt = -T^^Pxx + V{x, t)tlj, (2.2) 

where S{x,t) is a complex function, h is Planck's constant divided by 27r, m is the mass 
of the particle and V{x,t) is the potential energy function. The subscripts denote partial 
derivatives. The result is a quantum Hamilton- Jacobi equation ^, [lO, 11| 



where we recognize the classical Hamilton- Jacobi equation on the LHS. On the RHS is an 
additional non-classical term that can be referred to as a "quantum potential". Our aim 
is to solve eq. fi2.3p using the method of characteristics. This method is usually applied 
to the solution of first order PDEs. The underlying idea of the method of characteristics 
is to transform a single PDE (or a set of PDEs) to a set of ODEs that are solved along 
a characteristic curve. The characteristic is defined by setting a dependence between the 
independent variables of the PDE (or PDEs), in our case x and t. ZEVCA makes the 
simplest possible choice for the characteristics: 

^ = ^ x{t) = x(0), (2.4) 

i.e. the characteristics are trajectories that remain at a constant position. The time depen- 
dence of the phase along the trajectory is given by inserting S{x,t) in the Lagrangian time 
derivative defined as 

d _ d dx d , , 

di ^di^dtdx' ^ ■ 

The choice we made in eq.f l2.4p equates the Lagrangian time derivative and the partial time 

derivative ^ = ^- The result of the substitution is 

where we have used eq. (12.30 . The integration of eq. (l2.6p requires Sx[x{0),t] and Sxx:[x{0),t]. 
Note that the initial position x(0) acts as a parameter. Defining 



Sn[x{0),t] 



(2.7) 

x(0),t] 



dx^ 

we can write a general equation of motion for ^^[^(O),^] Vn in the following manner. We 
take the n*^ spatial partial derivative of eq. (12.31) 

iSt)n + 7^(^l)„ + K = 7^^n+2, (2.8) 

2m 2m 



and insert the result in the time derivative of Sr, 



(2.9) 



where we have assumed that the time and spatial derivatives are interchangeable. Equa- 
tion ( 12. 9p reveals that the equation of motion of S'„ depends on the result of two sub- 
sequent equations (by the dependence on 5'„+2) and on all prior equations by the term 
{Sf)n = X]j=o (j)'S'j+i'S'n-j+i- Hence, the general characteristic solution of eq.f l2.3p yields an 
infinite hierarchy of ODEs defined by eqs. (12.91) where n = 0, l,...,oo. In other words, we 
have converted eq. (12. 3p to an infinite set of local but coupled ODEs that generate the time 
dependence of the complex phase and its spatial derivatives at position x(0). An approxi- 
mation can now be obtained by truncating the set at some n = N; this is done by setting 
Sn+i = Sn+2 = 0. We summarize the equations of motion of the ZEVCA approximation 



dSn 

IT 

Sn+1 



7^5„+2 - 7^iS!)n - K[x(0)]; n = 0, 1, AT, 
2m 2m 



(2.10) 



N+2 



0, 



where we emphasize that the partial derivatives of the potential are taken at a fixed position 
x(0). The initial conditions of eqs. (l2.1UI) are given by 

9"ln[?/'(x,0)] 



dx"- 



(2.11) 



x{0) 



where we used the relation S{x,0) = —ih\rL[ip{x,0)] from ansatz (12. ip . After we solve set 
(I2.10p . the wavefunction at x(0) is given by 



■ilj[x{0),t] = exp 



-^o[x(0),t] 



(2.12) 



B. The ZEVCA imaginary time Schrodinger propagator 



The imaginary time Schrodinger propagator takes the form exp 



Ht 
2 



where 



H 



+ Vix). 



(2.13) 



2m dx"^ 

Applying this operator to an initial function il){x, 0) yields a "wavefunction" il){x, r) that is 
a solution of a Schrodinger-like equation 



2i!r{x, t) = ~l^'^xx{x, t) + V{x)'4){x, r). 



(2.14) 
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where r plays the role of time. The mapping of r to pure imaginary time i = — (hence 
the term — "imaginary time propagation") transforms the last equation to the form of the 
TDSE 

h'^ ~ ~ ~ ~ 

ihtpiix, t) = -^i'xxix, i) + V{x)ij{x, i), (2.15) 

where for simplicity we allow a slight misuse of notation tp{x, r) — > ip{x, t). At this stage we 
insert in eq. (l2.15p an ansatz identical to eq. fl2.ip 



ipi^x, t) = exp 



'^S{x,i) 



(2.16) 



and apply the ZEVCA formulation (section III Ap . This yields equations identical to 
eqs. (12.101) and (12. lip with the only difference being t — > t. 



III. MODEL APPLICATIONS 

A. Tunneling probabilities 

In this section we calculate tunneling probabilities form an Eckart barrier using a single 
trajectory propagation. The potential is given by 

= [cosHPxr ^'-'^ 

where D is the barrier height and 1/(3 gives an estimate of the barrier width. The initial 
wavefunction is a Gaussian wavepacket 



ip{x, 0) = exp 



,2 ^ ' ^ ^ 



-ao[x - Xcj + -Pc{x - Xc) + -7o 



(3.2) 



where 1 / ^/oo relates to the Gaussian width and 7o = — ^ In(^) takes care of the normaliza- 
tion. Xc and pc are the average position and momentum respectively. The initial conditions 
of eqs. (l2.10p are obtained by inserting eq. (l3.2p in eqs. (12. lip 

So{x,0) = iaoh{x - XcY + pc{x - Xc) + 'Jo, (3.3) 

5*1 (x,0) = 2iaoh{x — Xc) + Pc, 

S2{x,0) = 2iaoh, 

Sn{x,0) = 0, n > 3. 
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The final tunneling probability T for an initial wavefunction centered at Xc ^ and 
having pc> is given by 

T = lim T{t), (3.4) 



t— >oo 

where 



oo 



T{t) = / \^{x',t)\'dx', (3.5) 
Jo 

is the time-dependent tunneling probability. The integration initiates at a; = since this is 
the position of the maximum of the barrier. T(t) can be expressed by a time integration 
of the probability current at x = 0. We show this by first writing the quantum mechanical 
continuity equation 

pt{x,t) = -J^{x,t), (3.6) 

where p = l^/'p is the probability amplitude and the probability current J is given by 

j=Ac^(^t^^). (3.7) 
m 

By inserting p{x,t) is eq. (13.51) we can write 

POO 

T{t) = / p{x',t)dx' (3.8) 







oo pt 

dx' / pt{x',t')dt' 



Jo 

oo ft 







dx' / Jx{x',t')dt' 
Jo 



t poo ft 

dt' Jx{x\t')dx = I J{0,t')dt'. 
Jo Jo 



where in the third stage we used eq. fl3.6p and in the forth stage we changed the order of the 
integration and performed the spatial integration. Inserting ansatz (12. ip in eq. (l3.7p reveals 
that 



|2 



1 



J= IZ^mSi) = — exp 
m m 



-253(5, 







^ 3^(5i), (3.9) 

where we have used the notation defined by eq. (12.70 . Inserting eq. (l3.9p in the final result of 
eq.dM]) yields 

Tit) = l^\xp|^^^t|M)l|g^[5^(0,t')]dt'. (3.10) 

Note that 5*1 is readily obtained by the propagation of set (I2.10p . To calculate T{t) we need 
to set N and solve eqs. (l2.10p with initial conditions (12. lip at x(0) = 0. Inserting 5*0(0, t) 
and 5i(0,t) in eq. (l3.10p completes the derivation of the tunneling probability. 



We turn to the numerical results. The parameters we insert in eqs. (l3.ip and (13.20 are 
D = AO, (3 = 4.3228, ao = 307r, = -1.5, and Pc = V2mE where E = 20 and m = 30. All 
quantities here and henceforth are given in atomic units, hence ^ = 1. In fig {1] we illustrate 
the potential and the wavefunctions. Note that the initial wavefunction is located close to 
the barrier maximum, for reasons that we discuss below. 
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FIG. 1: Plot of an initial Gaussian wavefunction propagating in an Eckart barrier. At t = 1 where 
the wavefunction is clearly divided into a reflected part and a transmitted part. The depiction of 
the barrier's width is in proportion to that of the wavefunction, whereas the height of the barrier 
has no physical meaning. The arrow indicates the direction of the average average momentum. 
The parameters of the system are given in the text. 

In figures [2I^a) and[2]^b) we depict the approximations to |'?/'(0,t)p and T(t) for a series 
of values of the truncation order A^. Note that |?/'(0,t)p is equal to the exponential term 
in eq. fl3.10l) . The relative error between the exact tunneling probability and the asymptotic 
value of T{t) (see eq.(l33D) for A^ = 2, A^ = 6, and A^ = 10 is roughly 20%, 4% and 
0.5% respectively. Clearly, the numerical results converge quickly to the exact quantum 
mechanical result. In figures [3](a) and [3](b) we plot |'?/'(0,t)p and T(t) for a set of A^'s 
where we take Pc = 0- Although the relative error in the wavefunction ip{0,t) still converges 
uniformly, the relative error in T{t) does not: the errors are 1.5%, 6% and 0.8% for N = 2, 
N = 4, and A^ = 6 respectively. 
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FIG. 2: ZEVCA numerical results for |?/^(0,t)p ((a)) and the transmitted probability T(t) ((b)) 
for a series of values of the truncation order A^. The system corresponds to figlT] where the exact 
parameters are given in the text. The relative error between the exact tunneling probability and 
the asymptotic value of T{t) (see eq.(l33D) for N = 2, N = 6, and = 10 is roughly 20%, 4% and 
0.5% respectively. 

The ZEVCA formulation of tunneling probabilities has a significant restriction on its 
use. Since the ZEVCA trajectories remain fixed, the formulation is very sensitive to the 
initial conditions inserted in the equations of motion. The choice of x(0) must satisfy two 
properties: (1) The derivatives of the potential (l^[a;(0)], n = 1,...,N) must have at least 
one term that is significantly different from zero. (2) The initial wavefunction at x(0) 
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FIG. 3: Same as figlH expect that pc = 0. The relative errors between the the exact tunnehng 
probabihty and the asymptotic value of T{t) (see eq. (j3.4p ) for = 2, = 4, and N = 6 are 
roughly 1.5%, 6% and 0.8% respectively. 

('?/'[x(0), 0]) needs to be significantly different from zero. These restrictions ensure that 
the initial conditions "capture" the wavefunction and the potential's surroundings at x(0). 
These restrictions have prevented us from choosing an Xc in the negative asymptotic region 
of the Eckart barrier since such a choice would have produced ?/'(0, 0) —>■ 0. 
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B. Energy eigenvalues 



The ZEVCA imaginary time propagator allows for a simple calculation of the first energy 
eigenvalue of a bound potential. As in the pervious apphcation, the calculation requires just 
a single trajectory propagation. We start with a short derivation that demonstrates how the 
first eigenvalue may be calculated using imaginary time propagation. An arbitrary bound 
potential defines a set of eigenfunctions j = l,2...,oo that satisfy H(j)j{x) = Ejcpj, 

where Ej are the energy eigenvalues and H is given in eq. fl2.13p . The eigenfunctions can be 
used as a basis set for the expansion of an arbitrary wavefunction iIj{x) 



(3.11) 



where 



j=0 I 



1. If we define a time scale ti = ^ then the operation of the imaginary 



time propagator for r ^ ri on an initial wavepacket ip{x, 0) yields 



lim exp 

t2>ti 



Ht 



Hx, 0) 



lim exp 



Ht 



h 



^aj<pj{x) 



(3.12) 



lim > 

-r>ri ^-^ 

ai(j)i{x) exp 



aj(pj[x) exp 



h 



Ell 

h 



In the first stage we inserted eq. fl3.11l) . in the second stage we applied the imaginary time 
propagator on the eigenfunctions and in the last stage we applied the limit. Inserting the 
relation r = in the final result of eq. fl3.12l) and equating it with eq. fl2.16p yields 



hence, 



ai(f)i{x) exp 



S{x, t) 



2iEit 
2Eii 



exp 



ls{xrt) 

n 



ihln[ai(f)i{x)]. 



Taking the t partial derivative of the last equation yields 

Ei-- 



^Si{x,i), 



(3.13) 



(3.14) 



(3.15) 



where we recall that this relation holds for r ^ ri =^ \t\ ^ ^^i- As we mentioned in section 
III Bt eqs. (12.101) hold for the imaginary time propagator with the substitution oft ^ t. Hence, 



h ~ hdS h 



1 

2m ' 



V 



(3.16) 
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Since this equation holds for every choice of x(0) we need to propagate only a single trajectory 
for a sufficiently long (imaginary) time t to calculate Ei. 

Liu and Makri have used a Bohmian related formulation, the trajectory stability method 
(TSM), to calculate energy eigenvalues 12 1. TSM [3] emerges from conventional Bohmian 
mechanics by constructing a hierarchy of equations of motion for spatial derivatives of the 
phase and the amplitude of the wavefunction. In reference |12j], Liu and Makri use TSM 
for imaginary time propagation at constant-position characteristics, as we do here. But the 
modification of TSM for imaginary time propagation is non-unique; moreover, producing 
constant-position characteristics is quite an elaborate procedure that must be repeated at 
every time step. In contrast, the ZEVCA transformation from the Schrodinger real time 
propagator to the imaginary time propagator is unique and the fixed characteristics are 
obtained for all time simply by choosing ^ = 0. As we demonstrate, the energy eigenvalues 
obtained using ZEVCA are of the same accuracy as using TSM while the formulation is 
decidedly simpler. 

For the sake of comparison, we consider two of the potentials that were studied in reference 
12l |. The first is a quartic potential V{x) = + x^. The parameters of the Gaussian initial 
wavepacket (eq. fl3.2l) ) are = 0.5, Xc = I and Po = with m = 1. The fixed trajectory 
that we propagate is at x(0) = 0. The results are very robust with respect to the choice of 
the initial parameters. In figJUwe depict Ei (eq. fl3.16l) ) as function of r (r = ^i) for a series 
of values of the truncation order A^. The relative error between the exact energy eigenvalue 
{El = 0.8038) and the ZEVCA approximations are roughly 38%, 4% and 0.6% for A^ = 2, 
A^ = 4 and A^ = 8 respectively. For A^ = 16 the relative error reaches 0.1%. We see a clear 
convergence to the exact quantum mechanical result as a function of A^. 

The second example is a Morse potential V{x) = D[l — exp(— ax)]^. The parameters of 
the potential and the mass are D = 0.1745, a = 1.026 and m = 1837/2. These parameters 
correspond to the vibration of an H2 molecule. The parameters of the initial Gaussian 
wavepacket (eq.f l3.2l) ) are ao = 4.5924, Xc = 0.1 and po = 0. The fixed trajectory is 
positioned at x(0) = 0. As in the previous example, the results are robust with respect 
to the choice of the initial Gaussian parameters and the position of the fixed trajectory, 
x(0). In figl5]we depict Ei as function of r for a series of values of the truncation order A^. 
The relative error between the exact energy eigenvalue {Ei = 0.0098565) and the ZEVCA 
approximations are roughly 1.4%, 0.08% and 2 ■ 10~^% for A^ = 2, A^ = 4 and A^ = 6 
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Quartic potential 




FIG. 4: A comparison between the exact lowest energy eigenvalue of a quartic potential and 
the results obtained using a the ZEVCA imaginary time propagator in eq. (j3.16p with series of 
values of the truncation parameter A^. The potential function is V{x) = + x^, where the 
numerical parameters appear in the text. The relative error between the exact energy eigenvalue 
{El = 0.8038) and the ZEVCA approximations for A^ = 2, = 4 and = 8 is roughly 38%, 4% 
and 0.6% respectively. 

respectively. In this example we see a faster convergence to the exact result as a function 
of than in the quartic case. The reason is that the parameters of the Morse potential 
correspond to a small perturbation from a harmonic oscillator potential. For the harmonic 
oscillator it is possible to show that the ZEVCA approximation for A^ = 2 yields the exact 
quantum result. For anharmonic oscillators, the truncation at A^ = 2 incorporates only 
l^ [x(0), 0], j = 0, 1, 2 in eqs. fl2.10l) . which is equivalent to making a harmonic approximation 
to the potential (note that the results obtained for N = 2 for both the quartic potential and 
the Morse potential correspond to the energy eigenvalue of the harmonic approximation to 
the two potentials respectively). Since the parameters of the Morse potential correspond 
to a smaller perturbation from a harmonic oscillator potential than the parameters for the 
quartic potential, the convergence with A^ is faster. 
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FIG. 5: A comparison between the exact lowest energy eigenvalue of a Morse potential and the 
results obtained by using the ZEVC A imaginary time propagator in eq. (j3.16p with series of values 
of the truncation parameter N. The Morse potential parameters correspond to the vibration of an 
H2 molecule (the numerical parameters appear in the text). The plot is a zoom of a complete graph 
that initiates at r = 0, £^1 ~ 4.5 • 10^^. The relative error between the exact energy eigenvalue 
{El = 0.0098565) and the ZEVCA approximations is roughly 1.4%, 0.08% and 2-10-^% for = 2, 
A = 4 and A = 6 respectively. 

IV. SUMMARY 

We have presented ZEVCA, a new approximation for quantum dynamics calculations that 
is a cross between a grid method and a semiclassical method. The formulation was applied 
to the calculation of tunneling probabilities and low energy eigenvalues, with surprisingly 
good accuracy. The ZEVCA formulation has several advantages: (1) The derivation of the 
formulation is straightforward and the equations of motion are readily solvable by standard 
numerical software. (2) The ZEVCA approximation yields the solution of the TDSE at a 
fixed and predetermined position. This allows for easy application to quantum quantities 
that require only local information. (3) The ZEVCA formulation requires the calculation 
of the potential and its derivatives only at x(0) (see eqs. fl2.10l) ). This is in contrast to 
most semiclassical methods that propagate trajectories in configuration space and require 
the calculation of the potential (or its derivatives) at each time step. (4) No root search 
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is needed as in many other semiclassical methods. (5) The extension to imaginary time 
propagation is easily attainable. (6) By taking — > oo, ZEVCA formally gives the exact 
quantum result. Although we have not conducted a rigorous comparison of timings, for the 
applications presented here ZEVCA was found to be several orders of magnitude faster than 
the exact quantum calculation using a Split Operator method. 

Still, the numerical implementations reveals a number of limitations. First, the conver- 
gence to the exact result as a function of A^ seems to be asymptotic, in the sense that there 
is an optimal choice of A^. Second, the method has difficulty at nodal positions. Surpris- 
ingly, this is not a result of the original ansatz (eq.f l2.ip ) but is a limitation imposed by the 
condition of fixed trajectories. Since the trajectories are fixed, clearly they cannot cross. We 
have demonstrated elsewhere that making the ansatz (12.11) but allowing for contributions 



from crossing trajectories can produce interference effects such as nodes and oscillations [13|. 
Third, in section [III Al we have presented limitations on the relation between the position of 
the initial wavefunction, the derivatives of the potential and the choice of x(0). These limita- 
tions originate from the need to incorporate in the initial conditions of ZEVCA "sufficient" 
data for a successful propagation. 

An alternative to the procedure we have described is to discard ansatz (12.11) and construct 
a hierarchy of ODEs for the wavefunction iljn[x{0),t] itself, instead of for the complex phase, 
Sn[x{0),t]. We have explored this direction but found that it produces very poor results. 
In the case of an initial Gaussian wavepacket it is readily verified that the truncation for 
N >2 does not entail any approximation to the complex phase derivatives S'„[x(0),0] (see 
eqs. (l3.3p ). This is not the case for the derivatives of the initial wavefunction ?/'„[a;(0), 0] 
itself, which go to infinity as a function of A^ for every choice of x(0). This observation 
provides additional justification for making the replacement ip — > e*"^/^ in the first place. 

We are currently working on several extensions of the ZEVCA formulation. Specifically, 
we are exploiting the local properties of ZEVCA for the calculation of thermal rates 14 1 . 
Similar to the implementations of ZEVCA we have presented in this paper, the calculation 
of thermal rates requires only local information. This allows us to obtain thermal rates 
for one-dimensional systems using only two trajectory propagations. We aim at exploring 
both ZEVCA and its close relative, Bohmian Mechanics with Complex Action, in multi- 
dimensional quantum systems. 
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